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OPTIMUM STEP DESIGN FOR CENTERING OF PISTONS 
MOVING IN AN INCOMPRESSIBLE FLUID 
by Izhak Etsion* and Bernard J. Ham rock 
Lewis Research Center 

SUMMARY 

Hydrodynamic effects are analyzed for a stepped piston moving within a tight- 
clearance tube filled with an incompressible fluid. Together with the hydrostatic ef- 
fects that were analyzed in an earlier paper, a complete solution is obtained and an op- 
timum step design for centering of the piston is suggested. The axial speed resulting 
from an axial driving force is calculated, and some experimental results for pistons 
falling in a water-filled tube are presented. 


INTRODUCTION 

The operation of hydraulic equipment often involves the sliding of a piston within a 
cylindrical bore. Among the applications are hydraulic valves, pumps, and actuators 
(refs. 1 to 3); viscometers and timing devices (refs. 4 and 5); and hydrostatic extrud- 
ers (ref. 6). An important factor in all these applications is the rate of fluid leakage 
past the piston. It can be shown that the minimum leakage occurs when the piston is 
concentric within its conduit and that the leakage past the piston increases as the second 
power of the eccentricity. For a fully eccentric, plain cylindrical piston the leakage is 
2^ times that when the piston is concentric (ref. 4). Hence, to control the leakage, one 
has to control the eccentricity. For the leakage to be minimized, the eccentricity has 
to be minimized as well. 

Some methods have been suggested to reduce the eccentricity (refs. 7 to 9), but all 
these methods were analyzed or tested only with a stationary piston. This was based on 
the assumption that the hydrodynamic effects due to the axial motion are negligible as 
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compared with the hydrostatic effects due to the pressure difference across the piston. 

In reference 9 a stationary stepped piston is analyzed and an optimum step config- 
uration that minimizes the ratio of leakage to centering force is found. It is the objec- 
tive of this work to confirm experimentally the analysis of reference 9 and to find an- 
alytically the importance of hydrodynamic effects. 


hydrodynamic effects 

Step Pressure 

A stepped cylindrical piston in a cylindrical tube of the type being considered is 
shown in figure 1 . The tube is filled with a fluid, and the piston moves with a velocity 
U relative to the stationary tube. The leading edge of the piston has the smaller diam- 
eter so that the hydrodynamic pressure buildup around the piston produces a centering 
force (ref. 10), Based on this centering effect and the assumption that any other radial 
load acting on the piston is small enough to cause only small eccentricity, the problem 
can be treated as a one-dimensional case. Thus, the governing equation for a nontilted 
piston, where dh/dx = 0, is 



with the following boundary conditions (for the hydrodynamic effects); 


p = 0 at X = 0 and at x = L 
(All symbols are defined in appendix A. ) 

Since the film thickness is independent of x, the axial pressure variation must be 
linear, as shown in figure 2 where the hydrodynamic pressure on the step ('p ) is a 

V s7j^ 

function of 0 , To simplify the calculation, the coordinate system is fixed at the piston 
and the tube moves relative to the piston. 

From the continuity of flow at the step the following can be written; 

^ ^ , Jle (. 

2 12p Xg 2 12/i(L-Xg) 


This equation can be rewritten as 


2 



(Pb)^ / hj, h® . ,, 
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( 2 ) 


where 


h = C (1 + e cos 6) 
m m' ’ 




(3) 

(4) 


and 


Of = 


'M 


m 
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(5) 


(6) 


If the analysis is restricted to small eccentricity ratios e, the following approximate 
expressions can be written; 


hi = ci(l + 3e cos e) 
m m 


O' 


(7) 

( 8 ) 


Substituting equations (7) and (8) into equation (2) gives 

6mUL(q> - 1)0(1 - 0) 


(Pb) = 


h ^2 
'^m 


(1 - i3)(l + 3e cos 0) + /Sof^ (l + — cos 9 


a 
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where 


After rearranging terms, the step pressure (p ) becomes 

V S/j^ 


( 10 ) 


("«)h 


_ 6MUL(g - l)j3(l - P) , 


ci[l + /3(«3- 1)] \l + >^h 


cos 9} 
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where 


■V 


h 


3e 1 + /3(g^ - 1) 
1 + - 1 ) 


The location of the center of pressure can be obtained directly from the lo- 

11 

cation of the centroid of the typical pressure distribution shown in figure 2. Hence 
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“i 
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which yields the simple e^cpression 



(13) 


The location of the center of pressure is important for eliminating undesirable tilting 
moments. 


Centering Force 

The centering lorce (fig. 2) can be written as 



Substituting equation (11) into this equation gives 

W, = -6|uUr(-— \ ig ~ U/3(l - g) f’’ 

[l + -{\ I 1 + cos H 

In order that < 1, from equation (12) the following inequality must hold: 

3e [l + /3(a^ - li] < 1 + - 1) 

For Q! > 1, this is valid for any e < 1/3. Assuming yj^ < 1 and using reference 11 
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for the evaluation of the pressure integral give 


Wh 


= GiridVR 



(a - Dm - 13) 
jl + i3( 0!^ - l')| yjj 


(1 - rl) 


- 1/2 



(14) 


The centering force diminishes to zero as the piston approaches the concentric po- 
sition. Therefore, it is in that region of small eccentricity where one would like to 
maximize the centering force. If e is very small, it can be assumed that « 1 and 
the following approximation can be written; 


(1 - >h) 


■ 1/2 


1 1 2 

= 1 + - Vh 

2 


(15) 


By making use of equations (15) and (12), equation (14) can be written as 


W 


h 


= 077 pUR 





(16) 


where 


— (3{a - 


1)(1 - /3) [l + )3(o^ - 1)] 
[l + 0(0^ - l]^ 


(17) 


Flow Race 

The flow rate past the piston shown in figure 1 can be expressed as 



By making use of equation (7) and the tables of integrals found in reference 11, this 
equation can be integrated to give 

= tcURC^Qj^ (18) 


where 


5 



(19) 


% - 1 + 


l+^(a®-l)l' ^ ’'h 



2 

H (as was done in the previous section) it is assumed that y-^ « 1 and that equation 
(15) applies, equation (19) can be rewritten as 


Q, - a 


1 + Koi - 1 ) 
1 + - 1 ) 


( 20 ) 


Friction Force 


The friction force acting on the plug is 



The shear stress t in equation (21) can be written as 


(21) 


" h 2 Idxy 


( 22 ) 


Substituting equation (22) into (21) and integrating once over the length give the friction 
force as 



Substituting equation (11) and the film thickness equations (eqs. (3) and (4)) into this 
equation gives, after integration. 


p ^ 27TjuURL ~ 
'^m 


(23) 


where 


jS + 1 - ^ + 3(a - 1)^;3(1 - j3) 

|/„2 _ ^2 1 + /3(«2-1) 



(24) 
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which, for small eccentricity ratios e, becomes 


p ^ 1 + ;3(a- 1) ^ 3 (a- 1)^;3(1 - 0 ) 

" 1 + j8(q!^ - 1) 


( 25 ) 


The Optimum Step 

One of the main objects in piston design is to maintain the eccentricity as low as 
possible, thus minimizing the leakage past the piston and also diminishing the possibility 
of hydraulic lock (ref. 1). Minimum eccentricity can be achieved by selecting the par- 
ameters a and so as to maximize the nondimensional centering force Wj^. 

In figure 3 the nondimensional centering force Wj^ is plotted against the step length 
parameter /3 for various step heights. An interesting result is the existence of a max- 
imum for at Q! = 1. 6 and ^3 = 0. 26. This is different from the hydrostatic case of 
reference 9, where the centering force increases monotonically with the increase of a. 
However, the present optimum step configuration resembles the one described in ref- 
erence 12 for a flat, stepped slider, where the optimum values are a = 1, 87 and 
/3 = 0.28, _ 

Regarding the leakage, it is clear from equation (20) that the minimum for Qj^ oc- 
curs whenever a = 1 or jS = 1, which corresponds to a plain cylinder as expected. 
However, from equation (17) it is clear that for these values the centering force van- 
ishes. Thus, a plain cylindrical piston can assume any eccentricity up to e = 1 with 
the result of a significant increase in the leakage as compared with the concentric posi- 
tion (ref. 4). 

In reference 9 a design parameter for a stationary stepped piston is defined as the 
ratio between the nondimensional leakage and the centering force. An optimum step was 
found that minimized this parameter. The same procedure can be used for a moving 
piston. Thus, we shall look for a minimum of a parameter E given by 

E = 1 (26) 

\W4 ^{a - 1)(1 - p) 

The values of a and /3 that minimize the parameter E will serve as the optimum de- 
sign criteria. 

Figure 4 presents the variation of the design parameter E with respect to step 
height a and location j8. The most interesting result is that a large margin exists for 
an optimum design, ranging from an a of 1. 3 to 1. 7 and a /3 of 0. 3 to 0. 4; anywhere 
in this range, E is within 7 percent of the absolute minimum. This is the same range 
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as in the hydrostatic case of reference 9. Hence, the same step configuration is bene- 
ficial for both the hydrostatic and hydrodynamic effects. 


ELAPSED TIME OF MOVING PISTON 


For a piston moving at a constant velocity in a tube filled with fluid, the axial driv- 
ing force T is balanced by the friction force and the force due to the pressure differ- 
ence across the piston. Hence, 

T = ttR^p^ ^g + ^h 

where the hydrostatic friction force F^ is given in appendix B along with the other 
hydrostatic results. Substituting equations (23) and (BIO) into equation (27) gives 

T = 77R^Pg + JrRC^p ? + 277 mUL ^ Fj^ (28) 

^m 

If all the fluid remains in the tube, the axial velocity of the piston is given by 

ttR^U = Qg + Qh 

which, by substituting equations (18) and (B5), becomes 

ttRC^ - - 

(29) 

From equation (29) the pressure difference across the piston is 

- ^ s) <3«) 

m-'g 


_ 6pRL y /j 
\ 


and by substituting equation (30) into equation (28) the driving force becomes 


T = 277|uUL 







(31) 


The dimensionless parameters Qj^ and Fj^ as well as and F^ (see appendix B) 
are, for the range of optimum step design, of order of magnitude 1. Hence, the er- 
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ror involved in neglecting the hydrodynamic effects depends on C„/R. For most prac- 

_ 2 in 

tical applications, the value of C^^/R is less than 10 , thus justifying the neglect of 

hydrod 3 mamic effects. 

However, in some cases where the fluid is free to flow out of the cylinder, the 
pressure difference across the moving piston becomes very small, and the hydrostatic 
effects are less significant. Such a situation occurs for instance in gun recoil bearings 
(ref. 13), where in spite of low radial clearance the hydrodynamic effects are more im- 
portant than the hydrostatic effects. 

For the piston falling in a tube, as described in the following sections, all the 
terms containing C^/R in equation (31) can be neglected as compared to unity. Also 
the drivii^ force T is 

T = ffR^Lg(Pp - Pj) (32) 

Hence, equation (31) becomes 

g(Pn - Pf) = (33) 

Q 

m'^g 

When the velocity of the piston is constant as it falls down the tube, the time it takes the 
piston to fall a distance H is 


t 


H 

U 


which by equation (33) becomes 


6pRH 

g(p - Pf)Q 
m^'^^p r^g 


(34) 


As a check on the validity on this equation, if one assumes a plain piston instead of a 
stepped piston (impl 3 dng that in eq, (34) is unity), equation (34) is in complete 
agreement with the equation developed in reference 4. 


TEST APPARATUS AND PISTONS 

The test apparatus used in the experiments is shown in figure 5. The precision 
U-tube had a bore of 1, 5885±0, 003 centimeters with a wall thickness of 1. 12 centimeters 
and overall length of 120 centimeters. The tube material was Pso’ex glass that was op- 
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tically clear with graduations marked on it. Only one-half of the U-tube shown in fig- 
ure 5 was used in the tests. The tube was alined in the vertical position within 0. 25 
degree. 

Three of the pistons used in the tests are shown in figure 6. All the pistons were 
made of 17-4 PH stainless steel and hardened. The diameters of the pistons were ma- 
chined within ±0.0025 centimeter. Table I gives the dimensions of the pistons tested. 
The dimensions given in the table are shown physically in figure 1. Figure 6(a) cor- 
responds to piston 7, figure 6(b) corresponds to piston 2, and figure 6(c) corresponds to 
piston 3. 


TEST PROCEDURE 

Tap water was put into the U-tube, filling both ends. The top of the left half of the 
U-tube (fig. 5) was plugged. The first 13 cm of the tube shown in figure 5 was not used 
for measuring the time of the falling piston since acceleration of the piston occurs in 
this portion of the tube. Beyond the first 13 cm the piston falls at a constant velocity. 
The fall of the piston was timed with stopwatches capable of measuring to within 0. 2 sec- 
ond. Tests were repeated to assure the accuracy of the results. 


TEST RESULTS 

Table II shows the test results and how well these results agree with the theory. In 
this table a is defined by equation (5), ^ is defined by equation (10). and is de- 
fined by equation (B7). The theoretical elapsed time t it took the piston to fall 50. 8 
centimeters as described by equation (34) is also given in the table. The parameters in 
equation (34) and their respective values, which were constant for all the tests, are 

(1) Viscosity of fluid, \i, 0. 96x10 ^ N-sec/cm^ at 22° C 

(2) Radius of tube, R, 0.7943 centimeter 

(3) Distance piston falls, H, 50. 8 centimeters 

(4) Density of fluid, p^, 1 g/cm^ 

(5) Density of piston, p^, 7. 85 g/cm^ 

The parameter t” in table II is the experimental value of the elapsed time it took the 
piston to fall a distance of 50. 8 centimeters. 

The percentage difference between the experimental and theoretical values of the 
piston’s elapsed falling time can be written as 



In table II the value of A varies from -35 percent to +13 percent. 

In reference 4 it is noted that some of the plain cylindrical pistons wobbled as they 
fell down the tube. The same phenomenon was observed with piston 8 of the present set 
(an unstepped piston). The wobble was not apparent with unstepped piston 7, which had a 
clearance ratio C/R of 0. 0179 as compared with 0. 0264 for piston 8 . No wobbling 
occurred with any of the stepped pistons even at the large clearance ratios. This sug- 
gests a stabilizing effect of the step, as was observed in reference 4. For the stepped 
pistons, with the exception of piston 2, the agreement between theory and experiment 
was good, the deviation being 13 percent or less. The result for the nonwobbling, plain 
cylindrical piston also agrees with the theory. 


CONCLUDING REMARKS 

Hydrodynamic effects were analyzed for a stepped piston moving within a tight- 
clear ance tube filled with an incompressible fluid. Centering force, center of pressure, 
flow rate, and friction force were derived analytically. These results, together with 
the hydrostatic effects that were analyzed in an earlier paper constitute a complete sol- 
ution of the moving stepped piston and allow an optimum step design for centering of the 
piston to be calculated. The elapsed time of the moving piston resulting from an axial 
driving force was calculated. Experimental results for pistons falling in a water-filled 
tube were presented. The agreement between theory and experiments was good. For a 
plain piston with large clearance ratio (C/R > 0. 025), a wobbling motion was observed 
as the piston fell down the tube. The stepped pistons always moved smoothly without 
any wobble. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, September 16, 1976, 

505-04. 
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APPENDIX A 


A 



D 

D 

E 


1 

2 


e 


F 


F 


g 


F 


h 


g 

H 

h 

L 


P 

Pe 

Q 

Q 

Qh 

R 

T 

t 


t 


U 

W 


SYMBOLS 

100(t - t)/t 

larger radial clearance, (D - D2)/2 

smaller radial clearance, (D - Dj^)/2 

diameter of tube 

major diameter of piston 

minor diameter of piston 

design parameter, Q/W 

eccentricity 

friction force 

F /b-RC p 
g' m^e 

FhCm/27TiuURL 

gravitational acceleration 
distance piston falls 
film thickness 
length of piston 
pressure 

pressure difference across piston 
flow rate 

Qgl2M(LA')lC^p^ 

QhATOCm 
radius of tube 
axial driving force 

theoretical elapsed time of falling piston 
experimental elapsed time of falling piston 
velocity of moving piston 
centering force 

2Wg/37TLRPge 
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I 

i 

I 

I Wj^/97rMUR(L/C^)2e 

i 

X axial coordinate 

ot step height parameter, 

ji step length parameter, x /L 

O 

■ y defined by eq. (B2) 

& 

i defined by eq. (12) 

I e eccentricity ratio, 

^ 9 angular coordinate 

p viscosity of fluid 

Pj fluid density 

Pp piston density 

T shear stress 

Subscripts: 

cp center of pressure 

g hydrostatic 

h hydrod 3 mamic 

M larger clearance 

m smaller clearance 

s step 
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APPENDIX B 
HYDROSTATIC EFFECTS 


The hydrostatic step pressure as given in reference 9 is 

1 


(’s). = 


So 
a 13 


I + Pia - 1) 


1 + y COS 6 
g 


where 


1' = 3e 
g 


1-A 


1 + p{a - 1 ) 


O' - 1 


a 


The other parameters from reference 9 are 


r?) - 


+ i3 


Wg = i.LRPeWgC 


7tC 


Q = 

^ rr 


m 


PoQ 


12p 


e^g 


where 


^ a ;6(1 - /3)(a - 1) 
3 ^ 


[l + /3 (q'^ - 1) 


a 


Q = — 

g s 

1 + /3(a - 1) 

The hydrostatic friction force is found from 

'L /'Jr 


-// 


t„R d0 dx 


where 


(Bl) I 


(B2) 


(B3) 


(B4) 


(B5) 


(B6) 


(B7) 
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Integrating first over the length of the two portions of the piston gives 




gJ 


R d9 


which can be written as 


= 

g 


/' 


m 


(1 - o')(Pg) + ( 0 -+ e cos 0)p^ 


g 


R de 


Substituting (p_\ from equation (Bl) and integrating gives for small eccentricity 


F = 


0/^13 
1 + - 1 ) 


+ a 


(B8) 


(B9) 

ratios 

(BIO) 


where 
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TABLE I. - DIMENSIONS OF CYLINDRICAL PISTONS TESTED 


Piston 

Major diameter 
of piston, 

cm 

Minor diameter 
of piston, 

cm 

Length of 
major- diameter 
part of piston, 

^S’ 

cm 

Overall length 
of piston, 

L, 

cm 

1 

1.5606 

1.5545 

0.6299 

1. 5697 

2 

1. 5603 

1. 5494 

. 3150 

1. 5692 

3 

1. 5606 

1.5494 

.6299 

1.5710 

4 

1.5606 

1.5494 

. 9322 

1. 5697 

5 

1. 5606 

1. 5438 

.6299 

1.5687 

6 

1. 5469 

1. 5298 

. 6299 

1. 5690 

7 

1. 5608 


1. 5674 

1. 5674 

8 

1.5476 


1. 5672 

1. 5672 


TABLE II. - TEST RESULTS AND COMPARISON WITH THEORY 


Piston 

Step height 

Step length 

Dimensionless 

Smaller radial 

Theoretical 

Experimental 

Percentage difference 


parameter, 

parameter, 

hydrostatic 

clearance. 

time of 

time of 

between experimental 


a 


flow, 

Q 

g 

Cm- 

cm 

piston fall, 
t, 

sec 

piston fall, 
"t, 
sec 

and theoretical times, 
A 100 

1 

1.22 

0.4013 

1. 3652 

0.0140 

97.4 

102.0 

+ 5 

2 

1. 39 

.2007 

2. 0000 

.0141 

64.7 

42. 2 

-35 

3 

1.40 

.4010 

1. 6147 

.0140 

82.4 

77.4 

-6 

4 

1.40 

. 5939 

1. 3479 

.0140 

98.7 

104. 1 

+ 5 

5 

1.60 

.4016 

1. 8258 

.0141 

72. 9 

78.0 

+7 

6 

1.41 

.4015 

1. 6242 

.0208 

24.7 

27.8 

+ 13 

7 

1.00 

1.0 

1 

.0138 

136.7 

140.9 

+ 3 

8 

1.00 

1.0 

1 

.0204 

42. 4 

33. 1 

-22 
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STEP LENGTH PARAMETER, P 

Figure 3. - Centering force for various step configurations. 



Figure 4. - Design parameters for various step configurations. 




Figure?. - Test apparatus. 


(ai Piston 7 ismoothl. 


(bl Piston 2tXj/L = 0. 201. (cl Piston 3 (Xj/L = 0. ®l. 


Figure 6. - Examples of pistons used in experiments. 
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CONTRACTOR REPORTS: Scientific and 
technical information generated under a NASA 
contract or grant and considered an important 
contribution to existing knowledge. 


TECHNICAL TRANSLATIONS: Information 
published in a foreign language considered 
to merit NASA distribution in English. 

SPECIAL PUBLICATIONS: Information 
derived from or of value to NASA activities. 
Publications include final reports of major 
projects, monographs, data compilations, 
handbooks, sourcebooks, and special 
bibliographies. 

TECHNOLOGY UTILIZATION 
PUBLICATIONS: Information on technology 
used by NASA that may be of particular 
interest in commercial and other non-aerospace 
applications. Publications include Tech Briefs, 
Technology Utilization Reports and 
Technology Surveys. 


Details on the availability of these publications may be obtained from; 

SCIENTIFIC AND TECHNICAL INFORMATION OFFICE 
NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 

Washington, D.C. 20546 




